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ABSTRACT 


We calculate the one-loop QCD corrections to t > ix, using dimensional reduc- 
tion scheme, including QCD and supersymmetric QCD corrections. The analytic 
expressions for the corrections to the decay width are given, which can easily be 
extended to t > XT By. The numerical results show that the correction amounts to 
more than a 10% reduction in the partial width relative to the tree level result. We 
also compare the corrections in the no-mixing stop case with those in the mixing 
stop case. 
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1. Introduction 


The top quark has been dicovered by the CDF and DO Collaborations at the Fermilab 
Tevatron|l]. In the Standard Model t —^ W* + b is the dominante decay mode. Beyond 
the SM, in addition to the top decay into possible charged Higgs bosons plus bottom, a 
potentially important decay channel of the top quark is the supersymmetric decay into a 
lighter stop plus a neutralino, which has been extensively discussed at tree level|2]. It is 
generally expected that the lighter of the two stops is significantly lighter than the other 
squarks because the large top quark Yukawa coupling drives the diagonal stop masses to 
small values and enhances the off-diagonal mixing of left-handed and right-handed stops, so 
the present squark mass limits would not apply to the lighter stop. The best current lower 
bound on the stop mass is 55GeV and comes from LEP, operating at /s = 130—140 GeV[3]. 
The DO experiments at the Tevatron have excluded the existence of a stop lighter than 
100 GeV, albeit under certain assumptions|4]. Since the lightest neutralino is the lightest 
supersymmetric particle the decay t — 11X? could occur in a reasonably large volume of 
the parameter space with a sizeable branching ratio|2]. The one-loop radiative corrections 
to both t —^ W*b and t > Htb have been calculated [5,6] but the radiative corrections 
to t > ix? and t > Yo have so far not been calculated. In this paper we present the 
calculation of the one-loop O(a;) corrections to the top quark decay into the lightest stop 
plus a neutralino, including both QCD and supersymmetric QCD contributions. Our results 


can be generalized straightforwardly to the decay t —^ Xf bn, where b is a light sbottom. 


2. Tree-level 


In order to make this paper self-contained we first present the relevant interaction La- 
grangians of the Minimal Supersymmetric Standard Model (MSSM) and the tree-level decay 
rates for t > hx. The interactions of top and stop with neutralinos and gluinos are given 


by the Lagrangians|7] 


Lit, x9 = —VDt(Li; Pr, + Ri Pr)Xjti + h.c., (1) 


and 


Li = —g,T"t(a; — bis) Gati "hes 


where 
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Lı; —.4jc080— Csih0, La; = —A; sind — C; cos 0, (5) 
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Here Sw = sin 0w, Cw = cosOw, Pr R = i F ys), and Nj; are the elements of the 4 x 4 
matrix N defined in Ref.|7], which can be calculated numerically. T° = A*/2 are the Gell- 
Mann matrices and 0 is the mixing angle between left- and right-handed stops which are 


related to the mass eigenstates t; in Eqs. (1) and (2) by 


ti _ [{ cos@ sin tr 
i C e (10) 


This rotation matrix, Eq. (10), diagonalizes the stop mass matrix|8] 


—m,(A; + ui cot B) 


M2 -— M? +m? + 0.35 cos(28)M7 
t MZ. + mi + 0.16cos(28)M7 


in (Al 41605) SH) 


where M, M, are the soft SUSY-breaking mass terms for left- and right-handed stops, u 
is the supersymmetric Higgs mass parameter in the superpotential, A; is the trilinear soft 
SUSY-breaking parameter, and tan 8 = v5/v; is the ratio of the vacuum expectation values 


of the two Higgs doublets. 


The tree-level Feynman diagram for the decay t —^ tx? is shown in Fig.1(a), and the 


tree-level partial decay width is given by 
1 


Po = geo endis mid) [c Ly Qn ems — mi) 


+4Re(Li,Rij)mumgo| (12) 


where A(x, y, z) = (x — y — 2)? — Ayz. 


3. Virtual corrections 


Since the conventional dimensional regularization violates supersymmetry, in our calcu- 
lation we will use dimensional reduction technique[9], which preserves supersymmetry, for 
regularization of the ultraviolet divergences in the virtual loop corrections, although there 
is only a small difference between the both schemes to first order in the QCD and weak 
couplings. In fact, in dimensional reduction scheme, at the one-loop level the e-scalars con- 
vert the dimensional regularization result to the result which would be obtained by simply 
performing the numerator algebra in four dimensions|9]. To regulate the infrared divergences 
associated with soft and collinear gluon emission we will give the gluon a small finite mass 
A which is legitimate for our purposes since the non-Abelian nature of QCD does not show 
up in this order. We will also adopt the on-shell renormalization scheme[10] in which the 
coupling constant and the physical masses are chosen to be the renormalized parameters. 
The finite parts of the counterterms are then fixed by the renormalization conditions that 
the quark and the squark propagators have poles at their physical masses. For the QCD 
and SUSY-QCD corrections to the decay t —^ Xx), which we are considering, only the top 
quark mass and the stop mixing angle in the bare coupling need to be renormalized. By 


introducing appropriate counterterms the renormalized amplitude can be expressed as 


Mren = —iv2u(X9)(aPn + bPr)u(t) (13) 
with 
1 1 z 
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where AQP 


expressions for which will be given below. ôL 


* 


1j 


Y-QCD ; . ; ; 
2? and A7.2! —9CP are the vertex corrections from the irreducible vertex diagrams, 


; and 0 Rj; are the shifts from the bare couplings 


to renormalized couplings and , as mentioned above, can be found by renormalizing the top 


quark mass and the stop mixing angle: 


The counterterms and the renormalization constants in Eqs.(14)-(17) are defined by 


and 
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(22) 


Calculating the self-energy diagrams for the top quark in Figs. 1(b) and 2(a) we obtain 
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where the sum over i(— 1,2) is implied, and 
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and 
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Here, A = t — yp + log 4r with yg being the Euler constant and D = 4— 2e is the space-time 


€ 


dimension. The color factor Cr = 4/3 for SU(3) and u is the ’t Hooft mass parameter in 
the dimensional regularization scheme. Similarly, from Fig. 1(c), 2(b) and 2(d), one finds, 


for the stop, 


624, = TL |- Fin) i 2 (59 = 2m? (Gans) 4 Gia» 
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We have fixed the wave function renormalization constants and the top quark mass coun- 
terterm by the on mass-shell renormalization scheme condition. The mixing angle countert- 
erm is fixed by the requirement that 90 exactly cancel the remainder of the sum of all the 
ultraviolet(UV) divergent terms in the square of the renormalized ampltude, insuring the 
UV finiteness of the physical observables. From this requirement we found that the mixing 


angle counterterm simply is the negative of the counterterm 0Z;»; that is, 


This condition insures that all the ultraviolet divergences will cancel in the virtual corrections 
to the decay width, as will be seen below, and is in agreement with Ref.[11]. 


The calculations of the irreducible vertex corrections from Fig. 1(d) and 2(c) results in 


AC * * ^ 
= T {(Li;Pr + Ri; Pz)|A + 43 


5 


--m2(2€o + 2C1; — C15 + C31 — Cos) + m2 (2Co + 28 
TC + C3) = mio (2Co -+ 2C Pr Ci, = Co -+ Cə3)] 
+2(L1;Pr + Ri Pr)mimgo(Cu — C12) } (=p, ki, A, me, mg ), 
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respectively, where the sum over i(= 1,2) is implied. In Eqs.(35) and (36) 
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where 8;; = ajb; — biaj, and Co, Ci; are the three-point Feynman integrals given in the 


appendices of Ref. [12]. 
The virtual correction to the decay rate is then 
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where Og” and ae USY-9CP are the UV divergent parts of the QCD and SUSY-QCD vertex 


corrections, respectively. These are given by 


z aC 
gee — 8gusv QCD _ c ^, 


and 0, 04, S) and ie are defined to be 
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and 


J 
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(39) 
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We have checked analytically that all the ultraviolet divergences indeed cancel in the virtual 


corrections to the decay width, but the infrared divergent terms presist. 


4. Real corrections 


As is well known|13], to cancel the infrared divergences in the virtual corrections one 


needs to include real gluon emission, namely, t — tix 8g, as shown in Figs.1(e,f). As above, 


we will regulate the infrared divergences associated with the soft and collinear real gluon 


emission by the same finite small gluon mass A. In the calculation of the corrections due to 


real gluon emission to the partial width, it was necessary to perform the integration over the 


three-body phase space. After tedious but straightfoward calculations we obtained 


asC p 1 
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Here we adopt the notation of Ref.[14] where the definition of the functions I;, Lij (Mmi, mz, , mee) 
can be found. We also have checked numerically that the infrared divergences in óL pea} and 


OD wire do indeed cancel. 


5. Numerical results and discussions 


In the following we give the numerical results for t — ¢,¥9, where X? is the lightest 
neutralino. In our numerical calculation we fixed M = 200 GeV, u = —100 GeV and 
we used the relation M’ = 545M [7] to fix M'. For the parameters in stop sector we 
assumed M;, = Mz, and took the combination A; + cot to be one parameter. Note 
that (A, + «cot 3) = 0 corresponds to the case of no mixing in the stop mass matrix, 
Eq.(11). There are then three free parameters in the stop sector and we chose m;,,tan 8 
, and (A, + ucot 8) as the three independent parameters. Other input parameters are 
mz = 91.188GeV, em = 1/128.8, and Gp = 1.166372 x 10 ?(GeV) ?. The W mass was 


determined from [15] 


2 
1 
m3.(1— mw) m 45 
wi M! vV2Gr1-— Ar e 
where, for a heavy top, Ar is given by [16] 
N, 2 2 
Ara Nomi ua) 
16r sy mMiy 


Figure 3 shows the relative correction to the decay rate d/o, Fo being the tree-level 
rate, as a function of the lighter stop mass assuming mg = 500GeV and tan 8 = 11. The 
solid curve in Fig.3 corresponds to A; + cot 9 = 0, the no-mixing case, while the dotted 
curve corresponds to A; + u cot 8 = 100GeV, a mixing case. Note that in Fig.3 the lightest 


neutralino mass is ms 


x» = 68 GeV. It is clear that the correction in the mixing case is 


8 


larger than in the no-mixing case and can reach -20% for m;, = 100 GeV. Figure 4 shows 
the dependence of the relative correction to the decay width on the value of gluino mass 
for m;, = 50 GeV. Other parameter values are the same as in Fig.3. For the solid curve 
mj, = 50 GeV and m;, = 64 GeV and there are two peaks at m; = 112 GeV and m; = 126 
GeV due to the fact that m, — 176 GeV and the threshod for open top decay into gluino 
and stop is crossed in these regions. For the dashed curve m;, — 50 GeV and m;, — 194 
GeV and there is only one peak at m; — 126 GeV. When the gluino mass is heavier than 
200 GeV the correction in the mixing case is larger than in the no-mixing case and both 
corrections increase with gluino mass. Decoupling effects do not occur here, in contrast to 
the virtual SUSY corections to the decay and production processes in the SM. In Figure 5 
we present the dependence of the relative correction to the decay width on the value of tan 8 
assuming mg = 500GeV, m; = 50 GeV and A; + cot = 100 GeV. Only in the region 
where tan  « 2 is the correction to the decay width very sensitive to the value of tan f. 

In conclusion, we have shown that the one-loop QCD and SUSY-QCD corrections to 
i tx? can exceed -10% of the tree level partial width in both the no-mixing and the 


mixing case of stop masses, and these corrections are not sensitive to tan 9 for tan B > 2. 


Note added: While preparing this manuscript the preprint of A.Djouadi, W.Hollik and 
C.Junger (hep-ph/9605340) appeared where the QCD correction to the process t — hx? is 
also calculated. But Eq.(14) of their original paper were not correct. Very recently, in their 
revised version this mistake has been corrected by them. We thank W.Majerotto for useful 


communication. 
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Figure Captions 


Fig.1 Feynman diagrams for the tree-level process t — 14 x) and the QCD corrections. 

Fig.2 Feynman diagrams for the SUSY-QCD corrections. 

Fig.3 The relative correction ôI /To to the decay rate as a function of the lighter stop 
mass assuming mg = 500GeV and tang = 11. The solid and dotted curves correspond to 
A, + cot 3 = 0 (no mixing) and A; + ui cot 8 = 100GeV (mixing), respectively. 

Fig.4 The relative correction ôl /To to the decay rate as a function of the gluino mass 
assuming m; = 50GeV and tan = 11. The solid and dotted curves correspond to A; + 
ii cot 3 = 0 (no mixing) and A; + wcot 8 = 100GeV (mixing), respectively. 

Fig.5 The relative correction ôl /To to the decay rate as a function of tan 8 assuming 


mg = 500GeV, mz, = 50 GeV and A; + ui cot 8 = 100GeV. 
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